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Using microlensed quasars to probe the structure of the Milky Way 
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ABSTRACT 

This paper presents an investigation into the gravitational microlensing of quasars by stars and 
stellar remnants in the Milky Way. We present predictions for the all- sky microlensing optical 
depth, time-scale distributions and event rates for future large-area sky surveys. As expected, 
the total event rate increases rapidly with increasing magnitude limit, reflecting the fact that 
the number density of quasars is a steep function of magnitude. Surveys such as Pan-STARRS 
and LSST should be able to detect more than ten events per year, with typical event durations 
of around one month. Since microlensing of quasar sources suffers from fewer degeneracies 
than lensing of Milky Way sources, they could be used as a powerful tool for recovering the 
mass of the lensing object in a robust, often model-independent, manner. As a consequence, 
for a subset of these events it will be possible to directly 'weigh' the star (or stellar remnant) 
that is causing the lensing signal, either through higher order microlensing effects and/or 
high-precision astrometric observations of the lens star (using, for example, Gaia or SIM- 
lite). This means that such events could play a crucial role in stellar astronomy. Given the 
current operational timelines for Pan-STARRS and LSST, by the end of the decade they could 
potentially detect up to 100 events. Although this is still too few events to place detailed 
constraints on Galactic models, consistency checks can be carried out and such samples could 
lead to exciting and unexpected discoveries. 
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1 INTRODUCTION 

Gravitational microlensing is now established as a powerful tool 
for various aspects of Galactic science. The concept of gravitational 
microlensing, where a background object is temporarily brightened 
by the gravitational lensing effect of a passing foreground object, is 
believed to have been initially conceived by Einstein around 1912 
(Renn, Sauer & Stachel 1997). As he stated in a later publication 
(Einstein 1936]), his opinion was that "there is no great chance of 
observing this phenomenon". Fortunately the idea was later resur- 
rected by several authors (Liebes 1964, Refsdal 1964), in particu- 
lar |PaczySkr| ([1986]). Following the proposal of Paczynski, many 
experiments began to hunt for microlensing events towards the 
Galactic bulge and the Magellanic Clouds: O GLE ([Udalski et aT] 
|7992] >, EROS (|Aubourg et al.|1993| ), MACHO ( |Alcock et al.|1993] >7 
MOA (Bond et al. 2001 ). To date thousands of microlensing events 
have been detected toward several directions, including the Galac- 
tic bulge (Thomas et al. 2005 Sumi et al. 2006 ; Hamadache et al. 



2006 ) and even in surveys of bright nearby stars ( |Fukui et aJ7[ 2007 , 
Gaudi et al. 20081. 



2006 ), Magellanic Clouds (Tisserand et al. 2007 ; Wyrzykowski et 
al. |2009|[20T0l >, and M31 jCalchi Novati et al.|2005[|de Jong et al.) 
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The purpose of this paper is to deal with a previously under- 
exploited set of targets, namely quasars. Although quasar mi- 
crolensing has been observed for many years (e.g. |Chang & Refs^] 
|dal |1979||lrwin et al. |19891jAnguita et al.|2008) this has, with a 
few exceptions (e.g. |Byalko |1970| ), almost exclusively been lens- 
ing of background quasars by compact objects in the intervening 
(strong) lensing galaxy. In contrast, we are interested in the lensing 
of background quasars by stars in the Milky Way. This is a partic- 
ularly prescient subject because imminent large-scale optical sur- 
veys, such as Pan-STARRS ( |Kaiser et al. |2002) and LSST ( |Tyson| 
2002] |Ivezic et al. ||2008] >, will monitor the whole available sky 
every few nights down to magnitudes which are sufficient to moni- 
tor several million background quasars. If large enough samples of 
such microlensing events could be accrued then it will be possible 
to use them to constrain models of the Galaxy. 

One reason why quasars are potentially important targets is 
that they break a number of the degeneracies known to inhibit 
the analysis of microlensing events from sources within our own 
Galaxy. In general, for lensing of Galactic sources it is difficult 
to determine the properties of the lensing population because all 
of the information about the event is contained in one observable, 
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namely the event time-scale. As a consequence, there are a number 
of degeneracies which need to be broken in order to determine, for 
example, the lens mass or distance. For sources within the Milky 
Way, in general, we do not know the distance or transverse velocity 
of the star which is lensed. For quasars the interpretation is made 
significantly easier because we can assume that the source is at in- 
finity and that the transverse velocity is negligible. In this case the 
lens mass is simply a function of the (known) time-scale and the 
(unknown) lens distance and velocity, 
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where D L and \i t are the lens distance and proper motion, respec- 
tively, and t E is the event time- scale. This makes the interpreta- 
tion of quasar microlensing events significantly easier than that 
of Galactic sources and opens up the possibility of carrying out 
model-independent mass determinations using high precision as- 
trometry from, for example, ESA's Gaia mission (Perryman et al. 
2001). Even without a measurement of the lens distance/velocity, it 
could also be possible to use higher-order microlensing effects such 
as the microlensing parallax (Gould 1992) to break the degeneracy. 
Clearly quasar microlensing has the potential to become an impor- 
tant avenue for probing both stellar and Galactic astrophysics. 

In this paper we carry out Monte Carlo simulations to address 
this issue, presenting an all-sky prediction of the rate and optical 
depth of microlensed quasars. The paper is organised as follows. 
In Section [2] we describe the model and theory. In Section [3] we 
present the main results, dealing with the time- scale distribution 
( |3.1| ), the predictions for the optical depth and event rate ( |3.2[ ) and 
also the prospects for determining the lens masses \33\ . We present 
a summary and discussion in Section|4] 



2 THE MODEL 

We utilise a Monte Carlo simulation to give an estimate of the 
rate and optical depth for microlensed quasar events. The details 
of the simulation are described in the following subsections. Our 
microlensing model is constructed following Kiraga & Paczynski 

i 1994 J and is similar to a number of previous works (Wood & Mao 
2005| |Han |2008) . 



2.1 Stellar distribution 

In this work we want to obtain an all sky estimate of microlensed 
quasars by the Milky Way stars, which requires a model for the 
stellar mass distribution. To do this we follow the model of Binney 
|& Tremaine | ([2008 hereafter BT08). Unless otherwise stated all 
parameters are taken from Model I in table 2.3 of BT08 and the 
solar radius is taken to be 8 kpc. It is immediately clear that the 
stellar halo's contribution is negligibl^] Thus, only bulge and disc 
stars are considered in this paper. 

The bulge is treated as a bar with its longest axis inclined by 
about 20° to the line from Sun to the Galactic center and the stellar 
mass distribution in the bar is given by 



1 Even for lines-of-sight out of the plane, where the relative contribution 
of the halo will be at a maximum, the total mass is negligible compared to 
the disc. For example, if we adopt the model of BT08 for the disc, and the 
model of de Jong et al. (2010) for the stellar halo then the halo contributes 
less than 4 per cent of the total stellar mass for any given line-of-sight. 
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where s = yjx 2 + (y 2 + z 2 )lq 2 . The bulge is truncated at a radius 
r bar . The three dimension Cartesian coordinates (x, y, and z) are 
shown in Fig. [T] For the bar parameters we adopt a bar = 1 kpc, 
Qf bar = 1.8, q - 0.35, r bar = 3 kpc, andp barj0 = 1.22 M /pc 3 . We 
have set the normalisation of the bar mass so that the bar has the 
same mass as the original bulge (0.518 x 10 10 M o ). 

The disc is modeled as the sum of two double-exponentials, 
which correspond to the thin- and thick-disc components. The 
Galactic disc density distribution in cylindrical coordinates is 



Pdisc(R,Z) : 
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(3) 



where S disc = 1773.12 M /pc is the central surface density, R is 
the radial distance, R A ^ C - 2 kpc is the disc scale-length (both com- 
ponents are assumed to have the same scale-length), zo - 0.3 kpc, 
and z\ - 1.0 kpc are the scale-heights for the thin and thick com- 
ponents, respectively, and a = 0.9333 and a\ - 0.0667 are the 
respective normalisations. The total mass of the disc in this model 
is5.13xl0 10 M o . 

2.2 Lens Mass Function 

For the mass function we use the four part power-law distribution 
model of |Kroupa |p002] l, 



6(M) ee 
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{ 0.4687M" - 3 
0.0375M" 1 - 3 
0.0187M" 22 

10.0187AT 2 - 7 



0.01 < M/M < 0.08, 
0.08 < M/M < 0.5, 
0.5 <M/M < 1.0, 
1.0 < M/M < 120. 

The normalisation k is obtained through the integral, 

J r 120 
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(5) 



where Pdiscsoiar is the stellar density in the solar neighbourhood 
(which can be derived from equation [3}. We find that k - 
0.762 M" 1 pc -3 . Note that unlike some previous authors we do not 
adopt different mass functions for different populations; as argued 
by |Bastian, Covey & Meyer ] ( |2010| ), there is no clear evidence to 
support strong variations in the initial mass function. 

We invoke the following simple approach for dealing with 
stellar evolution. The objects with initial masses 0.01 ~ 0.08 M© 
and 0.08 ~ 1.0 M© are assumed to become brown dwarfs and 
main-sequence stars, respectively. The stars with initial mass 1.0 ~ 
8.0 M are assumed to evolve into 0.6 M© white dwarfs. The stars 
with 8.0 ~ 40.0 M© are assumed to evolve into 1.35 M© neutron 
stars, and more massive stars are assumed to evolve into 5.0 M© 
black holes. 

We assume that the lens number density is proportional to 
the stellar mass density given in equations ^ and (3]) (i.e. n oc 
Pbar + Pdisc)- For convenience, we rewrite the total stellar distribu- 
tion equation as follows, 



Ptotal - Pbar + Pdisc - Pdisc,solarX(^ ^» AJ- 



(6) 



From this we can deduce the lens number density within a mass 
range, M, M + dM and a velocity range v, v + dv, as follows, 



/!(/, b, D L , M, v u v b ) = xft b, D L ){ (M)f(v h v b ). 



(7) 



where /, b are the galactic coordinates and D L is the distance to 
the lens, £ is the lens mass function given in equation ^ and 
X(/, b, D L ) is the dimensionless mass distribution (see equation[6|. 
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2.3 Lens luminosity 

In addition to the lens mass, we also consider the lens luminosity. 
This allows us to investigate the possibility of detecting the lenses 
and measuring their proper motion and distance, which will break 
the degeneracy in equation {T} and allow a model-independent mea- 
surement of the lens mass. 

For the main- sequence stars we estimate the lens bright- 
ness using the mass-luminosity relation of |Cox | ( |1998| ), while 
we assume that all non-main-sequence lenses are dark. The mass- 
luminosity relation of Cox (1998) is given in the V-band in the 
Johnson-Cousins UBVRI system. In order to compare with LSST 
observations (see Section [2!6| ) we need to transform the V-band lu- 
minosity into a luminosity in the LSST ugrizy photometry system. 
To do this we use the empirical colour transformations given by 
|Jordi et al | ( [2006l ). 

As the lenses are distributed at different distances from the 
Sun, we need to carefully consider the Galactic dust extinction. To 
do this we utilise the Galactic dust map of Schle gel et al.| |l998 
hereafter SFD), combining this with a model of the interstellar 
medium from BT08, as follows. We first construct a model to cal- 
culate the total extinction along a given line-of- sight for the photo- 
metric band x, 



A x (oo) 



= k x f 

Jo 



p g dl, 



(8) 



where k x is a constant (i.e. this does not vary across the sky) and 
the distribution of the interstellar medium (p g ) is modelled using 
equation (2.211) of BT08. We estimate k x by carrying out a least- 
squared fit between A x (oo) and the observed value of the total ex- 
tinction from SFD using 10 6 lines-of- sight uniformly distributed 
across the sky. Once we have estimated k x we can calculate the 
extinction that a lens is subject to using the following formula, 

A x (D L ) = k x p g dl. (9) 

Jo 

2.4 Velocity Distribution 

To calculate the event rate we need to know the velocities of the 
observer, lenses and sources. For simplicity we neglect the motion 
of the Earth around the Sun, the so-called parallax effect (Gould 
1992), returning to this briefly in Section [33] For the calculation 
of microlensing event rate only the relative lens -source transverse 
velocity (i.e. the transverse velocity of the lens, in the lens-plane, 
relative to the observer-source line of sight), \ t , is needed, which 
can be written, 



v t ~ v L - v for D L <<c D s 



(10) 



where D L and D s are the distances to the lens and source, respec- 
tively, and v and v L are the components of the velocity tangential 
to the line-of-sight for the observer and lens, respectivel)]^] 

We model the velocity distributions of the lens stars in the 
Milky Way using Gaussian distributions, adopting the Solar mo- 
tion with respect to the local standard of rest from Schonrich , Bin-| 
|ney & Dehnen | ( |2010| ). We take two separate tri-variate Gaus- 
sians for the thin and thick discs, using the means and disper- 
sions from table 2 of |Smith et al. |{2007) . For the bar stars, the 

2 Note that although our sources are at cosmological distances, the trans- 
verse velocity equation reduces to the simple geometrical one given in 
Equation flO) for lenses residing at non-cosmological distances (see, for 
example, appendix B of Kayser, Refsdal & Stabell 119861. 
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random velocities are assumed to have Gaussi an distributions with 
(ay/, ay, ay/) = (113.6, 77.4, 66.3) km s" 1 ( |Han & Gould |1995) 
along the three axes of the bar. Here, the coordinates (jc", y", z") 
are centred at the Galactic centre; x" represents the longest axis of 
the bar, z" is directed towards the north Galactic pole. The bar ve- 
locity dispersion needs to be computed in the Galactocentric cylin- 
drical coordinate system by, 

en) 

(12) 

(13) 

where (o- x ,o- y ,o- z ) = (110.0, 82.5, 66.3) km s _1 are the velocity 
dispersions along the x, y, z axis as shown in Fig. [T] which are 
calculated from (ay/ , cry , ay/) given the assumed bar angle of 20°. 
The resultant velocity dispersion is a function of the azimuth qp. 
The rotation component of the bar velocities v^bar are estimated 
following [Han & Gould |(I995) , 

v _ hmaxij^) R < Ikpc, solid body rotation, 
qp ' * \ Vmax R > \ kpc, flat rotation. 

where v max = 100 km s _1 is adopted. 

For the calculation of the relative lens- source transverse ve- 
locity, \ t , we need to convert the lens velocity components from 
Galactocentric cylindrical coordinates into solar-centric spherical 
coordinates. The conversion is done as follows, 
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■ sin b , 



(15) 



where (v r , V/, v b ) are the velocity components in solar-centric spher- 
ical coordinates corresponding to the radial, galactic longitude and 
galactic latitude components, respectively. The angle a is the an- 
gle between the lines connecting the Sun, the projected position 
of the lens on the Galactic plane, and the Galactic centre (see 
Fig. [I}, which can be computed by, 



(16) 



a = sin -1 R sinl/R, R = ^R 2 + d 2 - 2R d cos /, 
where d = D L smb. 

2.5 Quasar Luminosity Function 

Quasars are very energetic and distant galaxies with an active galac- 
tic nucleus. By assuming an isotropic distribution, the distribution 
of quasars can be given by the so called quasar luminosity function 
(QLF). A lot of work has been carried out on the QLF in recent 
years ([Boyle et al. ||2000[ |Richards et al~||20061 |Fontanot et al 



2007 j Hennaw reTal. |2007URichards et al. |20091|Croom et al. 



2009 ), especially with the advent of surveys such as the Two De- 
gree Field (2dF) and Sloan Digital Sky Survey (SDSS). 

In this paper we assume that the quasars have a constant sur- 
face number density and follow the method of |Li et al~| ( |2007| ) to 
give an estimate of this number density. We utilise the A cold dark 
matter (A CDM) model in a flat universe (i.e. Cl m $ + Cl^o = 1) with 
^m,o = 03, Q. A ,o = 0.7 and h = 0.7, where h is the Hubble constant 
in units of 100 km s _1 Mpc -1 . We then simply integrate the QLF in 
the redshift interval z < 8 to derive the surface number density of 
quasars in the i- and z-bandqj We denote the total surface density 

3 Note that quasars with z > 8 will have a negligible contribution to the 
total surface number density (see fig. 1 oflLi et al. 120071. 
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Figure 1. Position of a lens star seen from the Sun, where (Dl, /, &) are the 
galactic coordinates of the star, (R, qp, z) are the Galactocentric cylindrical 
coordinates, and Rq is the Galactocentric distance of the Sun. The vectors 
e r , e/, and e^ are the unit vectors in the solar-centric spherical coordinates. 
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Figure 2. The surface number density of quasars as a function of limiting 
apparent magnitude in i- and z-bands. 



of quasars down to a particular limiting magnitude as N(< i) and 
N(< z) for the two photometric bands. Fig. [2] shows the resulting 
distributions. 

Since quasars are extragalactic sources, some of the light emit- 
ted from them will be absorbed by the dust in the Milky Way. We 
take this into account using the extinction map of SFD. The moti- 
vation for including the z-band in our analysis is to see whether the 
reduced extinction for this longer wavelength increases the quasar 
microlensing event rate compared to /-band. 



2.6 Optical Depth and Event Rate 

The microlensing optical depth is defined as the probability that a 
particular background source falls into the Einstein radius of any 
foreground lens star. As has been mentioned above in Section [23] 
we assume that the quasars are distributed uniformly across the sky. 
If we denote the surface number density as n s (l, b, D s ), the opti- 
cal depth along a given line-of-sight can be calculated following 
|Kiraga & Paczyfiski |(T994] ), 
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where G is gravitational constant, D = D L (D S - D L )/D S « D L for 
£) L «: D s , and M\ and M 2 are the mass limits of the lens stars. Note 
that the optical depth is independent of n s {l, b, D s ). 

The event rate is the number of microlensing events that oc- 
cur per unit time, which we again calculate following Kiraga & 
paczynski 1(1994) , 

4 Vg r°° r Ds 

T(l 9 b) = -^— dD s n s (l,b,D s )D 2 s dD L xd,b,D L )D 1/2 

c Jo Jo 

rM 2 nn 

dMto(M)M l/2 dvidv b v t f(v h v b ) 

JMi JJ 

■[ I dD s n s (l,b,D s )D 2 s r l 

Jo 

4VG r°s 



~M 2 



■ff 



dD L x(l,b,D L )D 1/2 dMto(M)M l/2 

J Mi 



dvidv b v t f(v h v b ). 



(18) 



The time-scale of a microlensing event, t E , which is defined to 
be the time for a source to cross the Einstein radius of the lens (e.g. 
|Mao 12008] ), is given by, 



^^V] 1 ' 2 for D L «D S . 



(19) 



where r E is the physical size of the Einstein radius in the lens plane 
and V* is the relative transverse velocity of the lens (equation [10|. 

During a time interval At, the total number of expected mi- 
crolensing events is trivially, 



AL 



: ^totaiA/ 1 , 



(20) 



dD L x(l,b,D s )D 



where r total is the total event rate integrated over all lines-of-sight 
and e is the detection efficiency. For simplicity, we assume e = 100 
per cent. This can be justified when one considers that the cadence 
of the surveys such as LSST will be of the order of a few days, 
which is suitable to detect practically all events in our time-scale 
range (see Section |3.1| >. We return to this issue in the following 
section. 



3 RESULTS 

3.1 Time-Scale distributions 

Using the model described above, we are able to investigate the 
time-scale distributions for quasar microlensing events. Fig. [3] 
shows the predicted event rate as a function of time- scale for the 
different Galactic components. The average time-scales for the 
bar, thin disc, thick disc and total (bar plus discs) are 31.7 days, 
44.4 days, 33.3 days and 41.2 days, respectively. The longest time- 
scales are obtained for thin-disc events towards / « 140 and / « 210 
degrees. In these regions the mean time-scale can reach as high as 
130 days; these events will be very well sampled and are likely 
to display microlensing parallax signatures (see Section [33) . Note 
that the average time- scale for disc-lensing events reaches a min- 
imum around the central parts of the Galaxy. The average in this 
region is reduced due to lenses on the far-side of the bar; these 
have large transverse velocities owing to their rotational velocity 
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acting in the opposite direction to the Sun's. The events with short- 
est time-scales are those caused by lenses in the Galactic bar, be- 
cause bar lenses have larger velocity dispersions and hence larger 
transverse velocities. Furthermore, unlike nearby disc-lenses, they 
do not share the disc's rotation velocity. 

The combined map in the lower-panel of Fig [3] closely re- 
sembles the time-scale map for the thin disc, since the majority 
of events are caused by thin-disc lenses (see Section [T2|. 



3.2 Optical Depth and Event Rate 

Arguably the most important quantity which we can calculate is the 
total event rate, i.e. the number of events which we expect future 
surveys will be able to detect. This is given in Table [T] for various 
apparent magnitude limits. Not surprisingly the total event rate in- 
creases rapidly with increasing depth, reflecting the fact that the 
number density of quasars rises steeply for fainter limiting magni- 
tudes (as shown in Fig. [2]). 

The contributions from the various Milky Way components 
are shown in Table [T] We can clearly see that the total event rate is 
dominated by the disc component, including the thin and thick disc 
components. The thick disc component contributes about 17 per 
cent of the whole disc events and therefore plays a non-negligible 
role in quasar microlensing. Although the total mass of the thin disc 
is 13 times more massive than the thick disc, the fact that the thick 
disc is less concentrated in regions of high extinction helps to boost 
its lensing signal. The bar component, with only around 10 per cent 
of the total mass of the disc, has significantly fewer events. As has 
been stated in Section |2.1[ the stellar halo will have a negligible 
contribution. 

The all-sky distributions of average optical depth and the event 
rates are presented in Fig. [4] The different rows of this figure show 
the average optical depth and the event rate for different appar- 
ent magnitude limits. As remarked above, the event rate increases 
with the increase of the apparent magnitude limit owing to the in- 
crease in source density. The event rate also increases with decreas- 
ing galactic latitude, tracing the scale-height of the discs. How- 
ever, as one reaches close to the Galactic plane extinction becomes 
important, effectively ruling out the ability to detect events with 
\b\ < 1.5°. Although one might expect this problem to be alleviated 
by operating at longer wavelengths, the gradient of the extinction 
(as a function of \b\) is so great that even in the z-band the event rate 
for the bar and disc only increases by around 30-50 per cent (see 
Table[l). 

The most crucial question which we can pose is how many 
quasar microlensing events do we expect future surveys to detect? 
There are a number of upcoming surveys which will map the entire 
visible sky every few nights. Two important ones are Pan-STARRS 
( |Kaiseretal.]2002l > and LSST { [Tyson |2002] ), both of which operate 
with a cadence of every few nights. We assume a 100 per cent de- 
tection efficiency, which should be a fair approximation since prac- 
tically all events will have durations longer than this cadence; from 
Fig. [3] we find that 95 per cent of all events have time-scales greater 
than 7.5 days (possible limitations to this assumption are discussed 
in Section]?]). Pan-STARRS will reach a /-band (z-band) magnitude 
limit of about iu m = 22.6 (zu m = 21.5) mag and will cover about 
30, 000 square degrees, while LSST will reach an /-band (z-band) 
magnitude limit of about ii im - 24.0 (zu m = 23.3) mag per visit 
and cover about 20, 000 square degrees. In the /-band we estimate 
that Pan-STARRS will detect three events per year, while LSST 
will detect around five events per year. Clearly these numbers de- 



pend on the detectability of events, which is an issue we will return 
to in the discussion. 

As mentioned above, conducting the survey in the z-band does 
not increase the number of detected events significantly. Further- 
more, for both of these surveys the z-band limiting magnitude is 
around one magnitude shallower than in the /-band, which cancels 
out any advantage gained from the lower extinction. As a conse- 
quence of these two factors, the total number of events which would 
be detected in this band is only 1 and 4 for Pan-STARRS and LSST, 
respectively, meaning that such a survey is optimally carried out in 
the /-band. 

In practice we believe that these are lower-limits to the num- 
ber of detected events, because the above calculation is based on 
events for which the lens and source are aligned to within one Ein- 
stein radius. This corresponds to a minimum amplification of 0.32 
magnitudes, which is easily within reach of these surveys with their 
milli-magnitude (or better) photometric precision. For example, if 
we allow for events within two Einstein radii (corresponding to a 
minimum amplification of around 0.06 mag) then the event rate will 
increase by a factor of 2. Although it will be challenging to identify 
events with such small amplification, it does means that surveys 
like LSST and Pan-STARRS could potentially observe around 16 
events per year. 

3.3 Lens mass determination 

It is important to investigate the characteristics of the lens popu- 
lation, particularly with respect to potential follow-up work on the 
lenses. 

In Section [2^2] we described our simplistic implementation of 
stellar evolution. From equation ^ it is trivial to obtain the rel- 
ative contribution of the different lens populations. This is sum- 
marised in Table [2] We see that by far the majority of lenses are 
main- sequence stars (~ 80 per cent), although a significant fraction 
are white/brown dwarfs. Around 1 per cent of events will be caused 
by neutron stars, which raises the remote, yet interesting, possibil- 
ity of quasar lensing from pulsars. As has been highlighted by pre- 
vious authors (e.g. |Dai, Xu & E samdin 2010), pulsar microlensing 
is of great interest because for these events it may be possible to 
determine both the distance and transverse velocity through radio 
observations (and hence, via equation]!] a model independent mass 
determination). 

However, the most promising avenue for carrying out model- 
independent mass determinations is for main-sequence stars. This 
is important because microlensing is the only method which can be 
applied to single stars (as opposed to all other methods of mass de- 
termination which are applicable only to stars in binary systems). 
As has been stated in the Introduction, to break the microlensing 
mass degeneracy requires a measurement of both the distance to 
the lens and its transverse velocity. However, to measure these lens 
properties requires the lens to be both bright and nearby. We in- 
vestigate this by plotting maps of the event rate for various lens 
properties in Figs.[5]to[8] 

In the near future instruments such as the ESO Gaia mission 
(Perryman et al. 2001 ) will measure parallaxes and proper motions 
to unprecedented accuracy for many millions of relatively bright 
stars. Its accuracy degrades rapidly for fainter stars, reaching ac- 
curacies of around 0.1-0.2 mas for stars at the magnitude limit of 
around G = 20 mag (B ailer- Jones ||2009| >. This corresponds to a 
parallax (and hence distance) error of around 20 per cent at a dis- 
tance of 2 kpc. Unfortunately, as can be seen in Fig. [5] we predict 
that very few lenses will be brighter than 20th magnitude in /-band. 
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Figure 3. The left-hand panels show the predicted microlensing event rate of quasars as a function of the time-scale. The right-hand panels show the corre- 
sponding maps of the predicted average event time-scale. We show the contributions of different Galactic components, the thin-disc, the thick-disc, the bar, 
and the combination of bar and discs from top to bottom, respectively. 
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rcyr- 1 ) 

mii m (mag) Bar Thin Disc Thick Disc Total 



21 

23 
24 
27 



estimates for /-band detection 

0.07 0.55 0.17 0.79 

0.51 3.53 0.86 4.90 

0.97 6.94 1.56 9.54 

5.13 34.75 6.92 46.80 



estimates for z-band detection 

21 0.12 0.91 0.23 1.26 

23 0.77 5.40 1.09 7.26 

24 1.50 10.31 1.95 13.76 
27 7.30 48.82 8.26 64.38 

Table 1. All-sky microlensing event rates for a range of magnitude limits. The different columns correspond to events due to the various Galactic components. 
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Figure 4. Distributions of the microlensing optical depths (top row) and event rates (bottom four rows). From left to right, the panels in each column 
represent maps of the bar component, the thin and thick disc components and the total contribution from all bar and disc components. Contours in the top row 
correspond to the log of the optical depth, while contours in the bottom four rows correspond to the log of events per year per square degree for the limiting 
source magnitudes shown on the right. 



In Fig. [6] we have plotted the map of event rate as a function of lens 
brightness and distance. From this we can estimate the total event 
rate for lenses bright enough and near enough to be within reach of 
Gaia. Equation {T} shows that the error on the lens mass will be en- 
tirely dominated by the uncertainty in the parallax; if we scale the 



Gaia errors linearly from the bright end to the faint end (i.e. from 
0.01 mas at G - 15 mag to 0.1 mas at G - 20 mag) and consider 
events for which we will be able to obtain a parallax error of 20 per 
cent, we find that we will obtain around 1 event every two years 
over the whole sky for quasars down to i - 24. If we allow for 
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Lens Type 


Event Fraction 




(per cent) 


Brown dwarf 


7.5 


Main sequence 


78.4 


White dwarf 


13.1 


Neutron star 


0.8 


Black hole 


0.2 



Table 2. The fraction of different lenses by population, following our sim- 
plistic implementation of stellar evolution (see Section [2!2) . These fractions 
are obtained trivially via the IMF (equation|4j. 



weakly magnified events (as above), this means that LSST will be 
able to detect around one event per year for which we will be able 
to robustly determine the distance (and hence lens mass) to within 
20 per cent accuracy using a parallax measurement from Gaia. 

Another potential source of accurate parallax measurements 
is the proposed SIM-lite space mission (Unwin et al. 2008 ), which, 
if it proceeds, will enable significantly better mass determinations. 
With an astrometric accuracy of 0.01 mas down to V = 20 mag, 
this will be able to follow-up a similar number of events to Gaia 
but with a ten-fold improvement in the accuracy. Even though the 
precision of the mass recovery will be greatly improved, the total 
number of events will not increase dramatically for SIM-lite since 
the limiting magnitude will be similar to Gaia; the only advantage 
to the event rate is the fact that the increased accuracy means it will 
be able to detect parallaxes for more distant lenses, which results in 
a doubling of the event rate compared to Gaia (i.e. around 2 events 
per year for LSST using i - 24 quasars). 

Even if we are unable to obtain accurate trigonometric paral- 
lax measurements, it will still be possible to accurately determine 
the proper motion for a significant number of lenses using ground- 
based astrometry. For example, LSST will be able to carry out as- 
trometry at around 3 mas for r - 24 mag. As can be seen from the 
lens luminosity function (Fig. [5} this covers a significantly larger 
number of events, many of which will have measurable proper mo- 
tion (see Fig. [7]). This would almost entirely break the lens mass 
degeneracy, but to obtain a unique mass measurement requires 
additional information, for example the distance (from photomet- 
ric/spectroscopic techniques) or through higher-order microlensing 
effects. One potential higher-order effect is microlensing parallax 
(Gould 1992), where the motion of the Earth as it orbits the Sun 
induces perturbations on the standard microlensing signal. Simula- 
tions of parallax frequency show that for distant sources the rate of 
parallax events towards the Galactic bulge can reach around 10 per 
cent (e.g. fig. 7 of |Smith et al. 2005 ). This fraction will probably 
increase for lines-of-sight away from the bulge, since the velocity 
dispersion of disc stars is smaller than bulge stars and hence the 
event time- scale will increase (and hence the parallax fraction will 
also increase as this is related to time-scale; e.g. see fig. 6 of Smith 
|et al. |2005] >. As is shown in Fig. [3] we find that 10 per cent of all 
events are predicted to have time- scales greater than three months. 
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Figure 5. The distribution of /-band lens magnitudes, after taking into ac- 
count extinction. Note that this corresponds only to the events caused by 
main-sequence lenses, i.e. a further 22 per cent of events have non-main- 
sequence lenses which we classify as 'dark' (see Table[2). 
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Figure 6. The all- sky event-rate map for luminous lenses as a function of 
lens distance and magnitude (after incorporating extinction). Contours cor- 
respond to events per year per magnitude per kpc for a limiting source mag- 
nitude of i = 24. Note that this corresponds only to the events caused by 
main-sequence lenses, i.e. a further 22 per cent of events have non-main- 
sequence lenses which we classify as 'dark' (see Table[2). 



4 CONCLUSION AND DISCUSSION 

In this paper we have carried out Monte Carlo simulations to make 
predictions for the microlensing of quasars by stars (and stellar 
remnants) in the Milky Way. Using a simple model of the Galaxy 
we have predicted the all- sky microlensing optical depth and have 
calculated the time- scale distributions and event rates for future 



large-area sky surveys, both for i- and z-band. We found that the 
total event rate increases rapidly with increasing magnitude limit, 
and predict that around 10 to 20 events will be detected per year 
from surveys like Pan-STARRS and LSST (see Section [T2| . Most 
events are caused by lenses in the thin disc, even though extinction 
hampers the detectability of low-latitude events. If we could iden- 
tify and monitor large numbers of quasars at low latitudes (e.g. [Lee] 
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Figure 7. Similar to Fig.[6] but for lens magnitude and proper motion 
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Figure 8. Similar to FigJSl but for lens distance and proper motion. Note 
that unlike Figs. [6] and |7p this plot also includes the non-main- sequence 
lenses (see Table[2). 



|et al. |2008| ), then it may even be possible to increase the number 
of detected events. 

We have also highlighted how such quasar microlensing can 
be used to carry out robust and model-independent lens mass de- 
terminations, predicting that with accurate astrometry from Gaia or 
SIM-lite it will be possible to accurately recover the lens mass for 
around one event per year. Microlensing is uniquely placed to pro- 
vide mass measurements for individual stars (as opposed to those 
in binary systems) and as such will undoubtedly play a crucial role 
in stellar astronomy. 

Although these results are very promising, there are a couple 



of technical challenges which must be overcome in order to detect 
these microlensing events. In this paper we have assumed a detec- 
tion efficiency of 100 per cent, which we have justified by noting 
that practically all events should be reasonably well-sampled by 
surveys like LSST. However, there are practical difficulties which 
will potentially reduce the efficiency. 

Firstly, it is well-known that quasars themselves are intrin- 
sically variable (e.g. Kozlowski et al. 2010). They can exhibit 
both long-term low-level variability as well as outburst phenom- 
ena. Indeed, one recent paper claimed to have detected a spectac- 
ular UV flare in a quasar behind M31 (Meusing er et al. ||2010| ), 
although the authors admit that an alternative explanation could be 
microlensing. Some authors have claimed that nearly all quasars are 
microlensed by intervening (cosmological) compact objects (e.g. 
Hawkins 1996), although more recent studies have challenged this 
hypothesis (e.g. |de Vries et al. |2005] >. This latter paper analyses 
a variability diagnostic (the structure function) of around 40,000 
quasars, concluding that the most-likely cause of quasar variability 
is flares due to disc instabilities. They remark that the variability 
is, in general, both chromatic and asymmetric. This is an impor- 
tant issue since the microlensing events which we are considering 
in this paper should be both achromatic and symmetric. Therefore 
one avenue for disentangling microlensing from intrinsic variability 
would be to utilise the multi-band photometry from surveys such as 
LSST. In addition, the amplitude of the quasar variability is gener- 
ally at a lower level than our microlensing amplitude (e.g. |Sesar et| 
|al. |2007| show that the rms scatter for quasars is ~ 0.07 mag in the 
g-band and, because the rms decreases for longer wavelengths, the 
variatbility will be even less for microlensing surveys carried out in 
i- or z-bands). 

Although most microlensing events are symmetric and achro- 
matic, this is not always the case. There are higher-order microlens- 
ing effects which can cause deviations from the standard symmetric 
and chromatic form, such as binarity in the lens or finite-source ef- 
fects. Both of these will probably be of only limited importance; for 
standard galactic microlensing the fraction of events which show 
clear binary signatures is of the order 3 per cent (e.g. Skowron et 
|al. |20 08, and references therein), while for finite-source effects it 
is even less. Although it is safe to assume that the fraction of bi- 
nary lens events for quasar sources will be similar to that for bulge 
sources, there is no reason to assume this will also hold for finite- 
source effects. If we assume that a typical quasar accretion disc is 
of the order 10 13 cm, then the angular size will be less than a few 
//as. This angular size is comparable to a giant star in the Galac- 
tic bulge, which means that the finite-source effect for quasar mi- 
crolensing should be no-more important than for standard Galactic 
microlensing; for these events finite-source effects only come into 
play when the amplification reaches around one hundred or more 
(e.g. |Cassan et al. |20 06). Therefore for the vast majority of quasar 
microlensing events the light-curves will indeed be symmetric and 
chromatic. For the limited number of events for which the magni- 
fication is sufficient to detect finite-source signatures, microlensing 
will enable constraints to be placed on the profiles of quasar ac- 
cretion discs, although this could only be done in very fortuitous 
circumstances (e.g. when a low-redshift quasar is very highly mag- 
nified). 

A second obstacle to detecting microlensing of quasars is the 
difficulty of identifying the quasars in the presence of a luminous 
lens. The quasars will be identified via cuts in colour-colour space, 
and so if the lens is sufficiently bright it could compromise the 
quasar selection. Fortunately in most cases the lenses will be sig- 
nificantly fainter than the sources (see Fig. [5}. We calculate that 
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for only a quarter of LSST events will the lens be brighter than 
the source and for just over half of all events the lenses will be 
fainter than the single-epoch magnitude limit of i - 24. However, 
this issue will be particularly important for the events from which 
we wish to determine the lens mass, since for these we require the 
lens to be luminous so that its distance and proper motion can be 
determined. For these events it will be important to spectroscopi- 
cally confirm the presence of a quasar behind the lens, or to obtain 
high-resolution imaging a number of years after the event so that 
the lens and quasar can be resolved. 

In contrast to the above challenges that hamper the detection 
of quasar microlensing events, there are avenues which could en- 
hance our ability to detect events. If we could build up an accu- 
rate astrometric map of the foreground lenses and their proper mo- 
tion (using, for example, data from the Gaia mission) then it could 
be possible to predict which quasars will be lensed in advance of 
the amplification. This method has been suggested as an approach 
to detect lensing by pulsars ( |Dai, Xu & Esamdin~||2010| > but it is 
equally applicable to any lens with an accurate proper motion mea- 
surment. Furthermore, this could be applied to fainter quasars, such 
as those provided by the complete co-added LSST survey, which 
will reach a depth of around i = 27. At this magnitude the number 
density of quasars climbs to around 3000 per square degree (see 
Fig. [2]) Another avenue is that of astrometric microlensing where, 
instead of detecting the photometric lensing signatures, we identify 
the astrometric shift of the image centroid as the quasar is lensed 
(e.g. |H0g, Novikov & Polnarev~|[l995l |Walker |[T^95] ). This will 
be a demanding task as it requires very high precision astrometry 
for our faint source population; the maximum astrometric shift is 
0.354 times the Einstein radius (Dominik & Sahu 2000 ) and for our 
events typical Einstein radii are around 1 mas. As above we could 
envisage targeting specific quasars based on predictions from lens 
proper motions, but this would still require astrometry at a level of 
at least tens of //as. 

In this work we have only considered lenses within our own 
Galaxy, yet one might wonder whether quasars could be lensed by 
stellar mass objects at cosmological distances. However, this is a 
remote possibility as it would require a very high density of such 
objects and even then the time-scales for such events are a couple 
of orders of magnitude longer than what we find for Galactic lenses 
(see, for example, Wambsgans s |20Q6) ]. The notable exception to 
this, as noted in the Introduction, is the phenomena of microlens- 
ing of macro-lensed (i.e. multiple-imaged) quasars by stellar mass 
lenses in the intervening lensing galaxy. Here the density of lenses 
is sufficient to produce an optical depth of order unity and hence 
detecting microlensing becomes feasible. Surveys such as LSST 
and Pan- STARRS will have an immense impact on this field as 
they will monitor the light-curves for many thousands of multiply- 
imaged quasars. 

Although we are unlikely to find events from lenses at cos- 
mological distances (with the exception of macro-lensed quasars), 
it may be possible to detect lenses from populations just outside 
our own Galaxy. For example, one could identify quasars behind 
the LM C (e.g. |KozIowski & KochaneF|2 009 ) or M31 ( |Huo et al. | 
2010) and study their lensing by foreground stars in these galaxies. 
As we mentioned above, there may have already been an event of 
this kind detected in a quasar behind M31 (Meusinger et al. |2010] >. 

In conclusion, we have shown that imminent ground-based 
surveys will be able to detect a number of quasar microlensing 
events every year. The efforts of Pan-STARRS and LSST will, after 
five years of operation, detect up to one hundred events. Although 
this is still too few events to constrain Galactic models due to the 



overwhelming numbers of parameters, it will allow for valuable 
consistency checks to be carried out on these models and could 
lead to exciting and unexpected discoveries. 
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